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We present the detailed formalism of the extremely correlated Fermi liquid theory, developed for 
treating the physics of the t-J model. We start from the exact Schwinger equation of motion 
for the Greens function for projected electrons, and develop a procedure by which this is expanded 
systematically in a parameter A, relating to the double occupancy The resulting Greens function has 
a canonical part arising from an effective Hamiltonian of the auxiliary electrons, and a caparison part, 
playing the role of a frequency dependent adaptive spectral weight. This adaptive weight balances 
the requirement at low uj, of the invariance of the Fermi volume, and at high uo of decaying as 
with a correlation depleted Co < 1. The effective Hamiltonian describing the auxiliary Fermions is 
given a natural interpretation with an effective interaction V e ff containing both the exchange Jij, 
and the hopping parameters tij. Simple but important shift invariances of the t-J model are noted 
with respect to translating its parameters uniformly. These play a crucial role in constraining the 
form of V e ff and also provide checks for further approximations. The auxiliary and physical Greens 
function satisfy two sum rules, and the Lagrange multipliers for these are identified. A complete 
and explicitly set of equations for the Greens functions to second order in A is given, satisfying 
various invariances. A systematic iterative procedure for higher order approximations is detailed. 
A superconducting instability of the theory is noted at the simplest level with a high transition 
temperature. 
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I. INTRODUCTION 



This work presents the detailed formalism of a newly developed framework for systematic calculation of the prop- 
erties of the t-J model, starting from the basic parameters t and J of the model. A subsequent paper Ref. (pQ) 
presents self consistent numerical results from the initial application of this theory, for the case of two dimensional 
square lattice relevant to cuprate superconductors. We will refer to extreme correlations as the limit U ^> t, so that 
the single occupancy constraint is enforced. The t-J model Eq. ( [Toj ) is the prime example of such a situation. In 
practice this theory applies already when U ^ Zt, where Z is the coordination number of the lattice. The theory and 
calculations presented are in the extremely correlated Fermi liquid (ECFL) phase discussed in Ref. (2). This phase is 
liquid like, and connects adiabatically to the Fermi liquid phase of weak coupling models such as the Hubbard model, 
while accommodating the extreme local interaction U —> oo. 

The t-J model has been the focus of intense effort for the last few decades, as reviewed in Ref. (3), particularly 
after its identification by Anderson in Ref. @ as governing many of the rich and complex set of phenomena High T c 
cuprates superconductors. The origin of the exchange part of the t-J model in an inverse expansion in the interaction 
U is familiar from superexchange theory. The relation J = is found there, so that large U leads to a small J. An 
early account of the model and the various sum rules can be found in the Ref. (5). More recently Zhang and Rice 
Ref. ([6]) gave an argument for reducing the three band copper oxygen model to an effective single band model, with 
a t-J form. Their method, apart from being more realistic, gives independent scales for t and J. 

Controlled calculations within this model are beset by two fundamental difficulties: (a) the non canonical nature 
of the single occupancy (Gutzwiller) projection of the electrons that changes the canonical anticommutation relations 
to a more nontrivial Hubbard (Lie) algebra and (b) the absence of any obvious small parameter for generating a 
systematic theory. The present author has recently formulated a method in Ref. ([2]) and Ref. (7) that overcomes 
these difficulties to a large extent. The basic idea is to approach the system starting from the limit of low particle 
density n = N e /N s (ratio of electron number to the number of sites), i.e. a generalized virial expansion. The density 
can be increased towards half filling systematically, as described below. Early applications to Angle Resolved Photo 
Emission (ARPES) experiments in Ref. (8) have been promising, and the general structure of the method already 
leads to non trivial and experimentally testable predictions in Ref. (9). This paper gives the details of the method 
and carries out a calculation to the lowest nontrivial order in a parameter A described below. The main ideas involved 
in this framework can be summarized as follows: 

• The Schwinger method: Ref. (|2J) utilizes the key observation that the Schwinger method dispenses with Wicks 
theorem, and replaces that step of canonical theory by a formal matrix (operator) inversion. The Schwinger 
equation for the Greens function typically involves a time derivative and a functional derivative with respect to 
a source potential V (defined more fully below). It has the great advantage over standard equations of motion 
in that the functional derivative generates all required higher order Greens functions. This is unlike say the 
BBGKY hierarchy, where one needs to import higher order correlations from elsewhere. For the t-J model 



Ref. ([2]) obtains a Schwinger equation described below in detail in in Eq. (40) and Eq. (41). Four our purpose, 
that equation may be illustrated schematically by the symbolic equation: 

(G -V)- y 1 - x)-g= (1- 7 ), (1) 
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where G 1 (Eq. (62 )) is a non interacting Greens function and 7 (Eq. (33)) is essentially the spatially localized but 
time dependent Greens function itself ~ Qiocai- Further Y\ is a Hartree type energy and X ~ (something) x 



contains the all important functional derivative with respect to V (both X, Y are defined in Eq. (41)). This is 
a convenient launching pad provided by Schwinger's method, since it is exact. However it is also intractable as 
it stands! There is no obvious small parameter, and the presence of the time dependent 7 on the right hand 
side represents the removal of states (and double occupancy) from the canonical theory and creates a new set 
of problems. We must understand and overcome these in order to create a practical and controlled scheme for 
calculations. We therefore push forward to the next set of steps. 

• Non canonical nature of the problem and its consequences: 

The non canonical nature of the problem is reflected in the 7 term on the right, which is a time dependent 
Greens function obtainable from Q itself (Eq. ([43])) . This 7 term contains an essential difficulty of the problem; 
it has a technical origin that we first discuss, and also an important physical aspect that we describe below. 

(*) Consider first the canonical theories, where one only has the 1 term in Eq. (HI). In order to get rid of 

— - - - - — _s_ 

sv 



the functional derivative operator X in favour of a (multiplicative) self energy, one uses X ~ (something) x 6 



to write: 

X g -^Sg = (something) x Q T Q, using —G = Q V Q, 

ov 

8 

following from T = — Q~ x , (2) 
ov 

wherein the vertex T is introduced. This gives the Schwinger Dyson relationship between the self energy E and 
vertex: 

E = (something) x Q T, so that (Gq x (m) - Yi - E) £ = 1. (3) 

This Schwinger Dyson construction necessarily requires that the vertex T reduce to unity at high frequencies, 
i.e. should be "asymptotically free". In case of the non canonical theory Eq. ([!]), a similar procedure fails. It is 
easily verified that the required good behaviour is lost because of the time dependent term 7 on the right hand 
side Ref. (10). The so defined vertex grows linearly with frequency, and invalidates the Dysonian self energy 
scheme. 

(*) The physical problem that is related to the non canonical 7 term has to do with the spectral weight of the 
projected electrons in a t-J model. Here basic sum rules give us insight into the origin, as well as a resolution 
of this fundamental problem. For non canonical electrons the high frequency behavior of the Greens function is 
^ 00 Juj w ^ n c o = 1 — f 7 rather than the familiar result for canonical electrons cq = 1. The depletion of cq from 
unity arises from the physics of single occupancy projection of the (non canonical) electrons q. Consider the 
relation Co = (c^cj^ -\-c\ a Ci a } , the process (c^cj^) suffers from the inhibiting requirement that the spin state a 
at site i must be unoccupied so that a double occupancy is not created by this process. On the other hand, if the 
numerator of Q(iuo) remains as Co at all frequencies, then the Fermi surface must enlarge in volume, and thereby 
violate the Luttinger Ward theorem of invariance of this volume 13 . We thus arrive at an appreciation of the 
tension between the conflicting requirements; at high frequency of fixing a known coefficient Co < 1, and at low 
frequency of a numerator almost unshifted from unity, for preserving the Fermi surface volume. A resolution is 
provided by the possibility of an adaptive (or smart) spectral weight, i.e the numerator of the Greens function. 
If a frequency dependent spectral weight can be found, so as to interpolate smoothly between the high and low 
frequency requirements, then both could be satisfied. 

(*)The above points suggest that the Greens function of the t-J model is usefully thought of as a 
product of two terms in frequency space i.e. 5^gx/i (Eq. ([42])) ? where g is a canonical Greens function and 
\i the caparison factor playing the role of an adaptive (or smart) spectral weight factor. The g term (i.e. the 
denominator) is required to be a canonical object with its poles and cuts as usual in a Fermi liquid, and defines 
the auxiliary Fermi liquid in this theory. The frequency dependent fi term (in the numerator) plays the role of 
the smart spectral weight; it reduces to the correct coefficient cq at high frequencies while recovering weight at 
lower frequencies. Thus a convolution in time domain into two suitable time dependent pieces could resolve this 



conundrum, and motivates the product ansatz in Eq. (42). 



(*)The [i term is also termed the caparison factor in Ref. (2), keeping in mind that it provides a second 
layer of dressing, over and above the dressing provided by the usual Fermi liquid type processes in g itself. This 
product ansatz is at the heart of the procedure described here and is seen to lead to a pair of exact equations 



for the two parts g and \i below in Eq. (65) and Eq. (66). 
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• Small parameter in theory: 

The t-J model is the sum of two highly non trivial terms, the kinetic energy projected to the space of single 
occupancy and the exchange energy It has no obvious small parameters making it especially difficult to deal 
with. Some inspiration is gained by examining the form of the analogous Schwinger equation for canonical 
theories, such as the Hubbard model. Again omitting details, the relevant equation can be written symbolically 
as: 

(GoHim)-UG- U-^)-G= 1, (4) 

where U is the Coulomb repulsion in the Hubbard model. Comparing with Eq. ([!]) suggests a simple approach 
to introduce a new parameter A. In its simplest form, we propose to study the modified problem symbolically 
expressed as: 

{Go\ft)-XY 1 -\X)-g.= (1-A 7 ), (5) 

with < A < 1, so that this equation Eq. ([5| interpolates smoothly between the Fermi gas and the t-J model. 
This appearance of the parameter parallels the way the Hubbard parameter U enters Eq. Q . The complication 
of the non canonical 7 term on the right is handled analogously to the Hartree term Y\. Unlike the repulsive 
Hubbard case, with an infinite interval [0, 00] for J7, the parameter A lives in a small and finite interval [0, 1]. 
The expectation is that low order perturbation expansion in A has a reasonable chance of capturing the physics 
of extreme correlations at A = 1. We show in the Appendix |A| that in the atomic limit, the role of A can be 
explicitly related to that of the fraction of double occupancy (and thus also density), so that tuning A smoothly 
adjusts this fraction between its two limits. 

• Effective Hamiltonian for the auxiliary Fermions with a pseudo potential: 

Keeping the caparison factor /i aside for a moment, we examine further the equations of motion for the auxiliary 
Fermion g following from Eq. ([5|, by an expansion in powers of A as in Sec. (II A). We would like to interpret 



these as the actual (canonical) equations of a suitable Fermi liquid, obtainable from a standard Hermitian 
Hamiltonian. However, we find that the equations do not immediately cooperate with this task. There are 
missing pieces that prevent such an interpretation at least immediately. The missing pieces must of course 
vanish identically at A = 1, since our starting point is exact in that limit. 

In order to resolve this issue, we construct an effective Hamiltonian H e ff for canonical electrons (fi a , fi a ), with 
the property that the (imaginary time) Heisenberg equation of motion for canonical electrons f ia = —[fi a , H e ff), 
match exactly the Heisenberg equation of motion for projected electrons Ci a = —[ci a: H t -j] : except for terms 
that vanish on enforcing the single occupancy constraint on the auxiliary f ia electrons. Thus we require 

[fia, H eff ] = ([c iCT ,# t -j])(£ ? £t)^(/jt) + (expressions involving /, /t that vanish at single occupancy). (6) 

We can then add these missing terms with (f,P) —> (c,$) to the Heisenberg equation of motion (EOM) for c 
and thereby obtain an auxiliary Fermi liquid that would be also "natural" , i.e. have all the standard properties 
of a Fermi liquicf^^l. One should therefore be able to use standard Feynman diagrams (Ref. (11)) to compute 
the properties of this auxiliary theory in powers of A, if one were so inclined. 



We find it straightforward to find such an effective Hamiltonian H e ff (Eq. (25 )) as described below in Sec. (II C ). 
A side remark on this construction is appropriate at this stage to give a physical meaning to H e ff. The 
kinetic energy of the projected electrons could also be written differently. An occasionally used alternate 
representation can be constructed by exploiting the fact that states satisfying the single occupancy constraint 
form an invariant subspace of a canonical Fermionic Hilbert space, under the action of a suitable set of operators. 
In this representation the Hilbert space continues to allow for double occupancy, i.e. is canonical, but the various 
operators representing the physical processes act only upon the singly occupied subspace, and produce states 
that are likewise singly occupied. Thus, with a = —a, we may write the kinetic energy part as: 

KE = ~ E t 1 " n -) flfi* (! " %-) (7) 

ij 

where n ia = fj a fi a - Since the exchange energy ^ - • J^ Si • Sj automatically conserves single occupancy, we will 
not write it out. The kinetic energy is thus a multi Fermi operator and represents both the propagation and 
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interaction between particles. To separate these functionalities, we introduce a parameter A here, it will turn 
out to be the same parameter as in Eq. (J5|, and write 

KE{\) = - £ tij (1 - A n ra ) flf ja (1 - A n js ) 

= " E *« fUia + A E *« ("» + n i*) + A " ^ ( 8 ) 

ij ij 

H d = ~^2Uj fLfja K^ja) , #d ^ 0. (9) 
ij 

The term if^ acts on the doubly occupied subspace and is null in the singly occupied space, and hence it may 
be dropped altogether. The remaining part of the kinetic energy term KE(X) has the structure of a four Fermi 



interaction between the canonical Fermions, and turns out to be a large part of H e ff in Eq. (25). 
In this representation, the physical electron Greens function Qij corresponds to the correlator —(((1 — 
n iai)fiai,fj aj 0- ~ n ivj)))i while — ((fi ai , f] aj )) wom d represent the auxiliary Greens function g(i, j). The capari- 
son factor \i seems hard to interpret in this language though. Also controlled calculations within this framework 
do not seem readily feasible. The ECFL formalism developed here presents a way to precisely glue together 
g and fx to give the physical Q. Its somewhat formal structure derives considerable benefit from the above 
connection; in particular Eq. (|8| helps in developing some intuition for g. For instance a physical interpretation 
of the auxiliary Fermions is provided by the fi a themselves, and thereby requiring the same number of auxiliary 
Fermions as the physical ones, as done below, is perfectly natural. 

(*) Invariances of H e ff and the emergence of the second Lagrange multiplier uq: 



In H e ff (Eq. (25)), the hopping parameter Uj is elevated to the role of an interaction coupling, in addition to its 
role a band hopping parameter. This feature needs attention, since we know that a constant (k independent) shift 
of the band energies — >• Sk +^t, or adding an onsite interaction through J^ — » J^ -\-5%j uj, is inconsequential 
for the t-J model, but it is not so for Eq. (8|. This "pure" gauge invariance is of primary importance in this 
kind of a theory, and must be addressed at the very outset to obtain a consistent and meaningful description of 
the t-J model. Such shifts could potentially lead to a change of the interaction strengths in H e ff, unless these 
shifts can be eliminated here too. This issue is addressed by first listing these shift symmetries of the model in 
Sec. pi]), and then requiring the approximation scheme to be shift invariant, at each order of A. 

Imposing the shift symmetries on H e ff has other important and interesting consequences. We shall see below 
that H e ff Eq. ( [25] ) is required to have a term with a Hubbard Coulomb like interaction, with coefficient u$. It 
ensures that arbitrary shifts of t and J do not change the physics of the problem, and can be absorbed into 
this uq parameter. Its magnitude remains undetermined by the symmetry operations. Along with the usual 



chemical potential /L£, it is fixed later by the two sum rules Eq. (87) and Eq. (88), one for the number of physical 



electrons and the other for the (identical) number of auxiliary electrons. 

In this work, we obtain a set of equations for the Greens function. These are essentially of the same form as in our 
recent earlier Letter Ref. (|2j), but differ in a few details due to the usage of the idea of the effective Hamiltonian and its 
shift invariances. The complete formalism contains two sum rules as already reported in Ref. ([2]), and now we exhibit 
two independent Lagrange multipliers available to enforce these. An iterative framework is carefully established, and 
calculations to second order in A are carried out. Higher order calculations should be feasible along these lines. To 
make the paper reasonably self contained, we give some details in brief that are already contained in our published 
work. 

The outline of the paper is as follows. In Sec. (II), we list the shift symmetries of the t-J model and obtain the 
exact equation satisfied by the Greens function. We also determine the form of the effective Hamiltonian H e f f for the 
auxiliary Fermions, such that the Heisenberg equations for the field operators are satisfied in a Hermit ian framework. 
In Sec. (III-IV) , we use the product ansatz for the Greens function to introduce and find the exact equations for the 
auxiliary Fermions and the caparison factor \i. In Sec. (V) we turn off the time dependent sources and write the 
exact momentum space relations between the self energy, the caparison factor and the physical Greens functions- 
these are the analogs of the Schwinger-Dyson equations for this problem. Appendix. (A) gives a detailed calculation 
in the atomic limit, and is expected to be useful in providing the reader some insight into the interpretation of the A 
expansion in terms of the number of doubly occupied sites. Section. (VI) summarizes the necessary equations needed 
for the next step discussed in Sec. (VII) namely the A expansion of various objects and the precise nature of the 
iterative expansion- which is the analog of the skeleton graph expansion in this theory. Several detailed calculations 
are needed to obtain the second order equations, and are detailed in Appendix. (B). Sec. (VIII) details the Ward 
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identities of this theory, which splits into two parts following the splitting of the Greens functions. Sec. (IX) gives the 
set of vertices defining the random phase approximation for this theory and Sec. (X) gives the formal results for the 
charge and spin susceptibilities within RPA and its low order expansion. Sec. (XII) concludes with some comments 
including a calculation of the superconducting transition temperature in this theory. 



II. THE t-J MODEL AND ITS SHIFT INVARIANCE 



We write the projected Fermi operators in terms of the Hubbard X operators as usual c ia —> X® a , c\ a — > X?° and 
c\ a ,Ci a —> Xf a . We study the t-J model given by 



= - £ t^xfxf -/x£xr + ^E J v { x f x T - xfxf) (io) 

i,j,o i,o ij,<7 

We will treat the two terms on an equal footing as far as possible. The statement of the model is invariant under a 
particular "pure gauge" transformation that we next discuss. Let us note the shift invariance of the two parameters 
in H. Consider the uniform (i.e. space independent) shifts of the basic parameters: 

t%j y tij Uf Sij , Jij y Jij -\- uj Sij , (11) 

with independent parameters u t ,uj. Under this transformation the Hamiltonian shifts as 

H^H+ (ut + ^uj^J N (12) 

where TV = ^2 icr X[ a is the number operator for the electrons. Let us note two simple theorems encoding this 
invariance: 

• Shift theorem- (I): A shift of either tor J can be absorbed into suitable parameters, leaving the physics unchanged. 

• Shift theorem- (II): The two shifts of t and J cancel each other when uj = —4 x u t . 



The first theorem is illustrated, for example, in the initial Hamiltonian Eq. (10) where it can be absorbed in the 
chemical potential /j, — > fi + u t , and serves to identify a second generalized chemical potential encountered later 
while making approximations. The second theorem is subtle as it leaves the chemical potential i± unchanged. It 
provides a measure of the equal handed treatment of t and J. We will find these almost trivial theorems of great use 
in devising and validating various approximation schemes later. 
In further work we need to add a source term via the operator A 

A=[ A(r)dr= £ f drVj^(T)Xj^(r)+ £ f dr ^{r)X^\r)Xf 2 (r), (13) 
Jo Jo -Jo 

with the usual imaginary time Heisenberg picture r dependence of the operators Q(r) = e rH Qe~ rH , and the Bosonic 
sources, Vj 1(T2 (r) at every site and also V^- 1<72 (r) for every pair of sites, as arbitrary functions of time. For any variable 
we define a modified expectation 

«Q(ti,t 2) ..)» = Tr [e-WT(e-*)} ' ( } 

with a compact notation that includes the time ordering and the exponential factor automatically. With the abbre- 
viation i = (Ri,Ti) for spatial R4 and imaginary time (r) coordinates, the physical electron is described by a Greens 
function: 

g„ i „ f [iJ} = -(W i Xf f0 ))- (15) 

Prom this, the variation can be found from functional differentiation as 

S 



sv * l( r Hn) ((Q(T2))) = «Q(r 2 )» ((Xj^in))) - ((Xj^irMr,))). (16) 

We note the fundamental anticommutator between the destruction and creation operators: 

{X^,Xp } = % («W 2 - (a lt r 2 ) X°^) . (17) 
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A. The Heisenberg Equation of Motion 



Let us now study the time evolution of the destruction operator through its important commutator: 



[X^,H] = -$>i [<W, - (aw) Xr*] + X°°* -nX°°' -Wjijia^) X°^ . (18) 



Here Jo = Jo is the zero wave vector, i.e. onsite exchange constant introduced for the purpose of satisfying the Shift 
Theorems (I- II) as explained next. The term in underbraces here and in the next equation, is introduced to make sure 
that the commutator reproduces the term with J^- — >• Jjj + uj 5{j correctly. Its necessity is seen when we note that 
under the transformation Eq. (12), the last term in Eq. (18) adds nothing, in view of the ordering of the operators as 



written, while the term with underbraces provides the correct transformation factor. Let us call this commutator as: 



[xr%H\ 



1 



Jo -fiX 



0(T l 



We next express the EOM for the Greens function in terms of A. 



(19) 
(20) 



B. Equation of motion for Q 

Let us compute the time derivative of the G. For this we need the derivative 

d n T{e- A X^{n)) = -T(e- A [Xt'(n),H])+T(e- A [A(r),Xt'}) 

[A(r),x^] = vr 2 (n)[x^(n),x^}-^2Kr 2 (n){xr°,x^} x^ 

3 



(21) 



3 3 

This follows from the definition of the time ordering and the form of A. Using this we find: 

dnSwAhft = - T f )5 hf (((s^ f - OiOfXr'))) + (([X^(n),H] X a /°(r f ))) 

-vr 2 {ri)Q^A^f]-Y.W 2 ^^f]-Y. v T 2 (w) {(xrHn)X^(n)X°/°(T f ))). 

3 3 

(22) 

We further use the abbreviations, 



(23) 



In terms of these, and using Eq. (19) we find the equation of motion in terms of A: 



(d Ti - riGwtii, f] = S[i, f]((S ai ,a f - v i o s X u i iUS )) 



+t[i,j] G^ f \j,f] - \jo Ga^lhfl + ((A ai (n)X a /°(r f ))} - vr 3 {n) 
-J2 v ?f 2g ^f&ft-I2 v ?r ^) {{X^{n)X°p{n)xy\r f ))). 



(24) 



We recall from the introduction, the discussion regarding suitably generalizing A of Eq. (20), in order to make 
connection with a Hermitian H e ff, and therefore turn to this task next. 
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C. Effective Hamiltonian 



We now construct an effective Hamiltonian of canonical Fermions that will turn out to govern the auxiliary Fermi 
liquid theory. The motivation for this construction is to cast the auxiliary Fermionic part of the ECFL theory into 
a natural and canonical framework, so that the equation for the g, i.e. the auxiliary piece of the full Q is obtainable 
from a Hamiltonian that is Hermitian and respects the usual Fermi symmetry of interactions under exchange. 

After some inspections a suitable Hamiltonian is provided by the expression: 

H eff = -£wU> + £(^o-M)/itJW + ^//, 

ij i 

Ve ff = \ S*«( ai<72 ) [{fljU + fj<nfj?i) ho J jo* + (h.c)] -\Y1 J v(w) fljj^fjoji^ 
ij ij 

+^«oM)44/ l52 / iff2 . (25) 

i 

with a Hermitian effective potential V^f = V e ff and assume no constraint on double occupancy for these auxiliary 
(canonical) Fermions fi a . The t and J parts reproduce the exact equations of motion as shown below with certain 
additional terms that vanish under the constraint of single occupancy. The parameter uq represent an effective 
Hubbard type interaction for these Fermions, giving a contribution uo^2i fj^fitfiifii- ^ s ma g n itude is arbitrary at 
the moment, since it disappears under exclusion of double occupancy. However it enables us to enforce the invariance 
in Shift-theorem- (I), where the shift of t and J can be absorbed in uq. It will turn out to play the role of a second 



chemical potential or Lagrange multiplier, in fixing the second sum rule Eq. (88). To illustrate this remark, note that 



adding a constant to t or J as in Eq. (12), adds an onsite four Fermi interaction term. In order to satisfy the Shift 
theorem -(I), we must compensate for this suitably, leading to the extra onsite term with coefficient ^0, which can 
absorb this shift. It is also verified that the Shift theorem- (II) is satisfied without the uo term. We emphasize that 
the uq term is both natural and essential for the purpose of satisfying the Shift theorem (I). Since the structure of the 
uo term is almost identical to that of we will most often "hide it" inside except when we explicitly display it. 
Thus unless explicitly displayed, we should read Jij — >• Jij — u Q Sij below. For analogous terms involving the X[ a 
operators as in Eq. (20), we can include uq in without any errors, since the uq term always vanishes due to the 



properties of these operators. 




FIG. 1: The pseudo-potential V e ff in the real space representation, where the wavy line represents Uj and the coiled line 
represents Jij. The first two interaction vertices have two undisplayed symmetric partners with the exchange i <H> j. 



Defining symmetric Cooper pair singlet operators 

VHi,i) = Y,°fU* = *fyi (26) 

with V^(j,i) = V^(i,j) we write 

V 'ff = lE** [(P f (M) +7^.7)) Hhi) + (h.c.)] -\Y, J ^ v\j ?V ( 27 ) 



ij 13 



9 



In momentum representation this reads 



ka 



4 N< 



(^1^2) W e ff(p 1 ,p 2 ]Ps,P4) fp^jLaJps^fp 



W e ff(pi,P2',P3,P4) 



°Pl+P2,P3+P4 



(28) 



where the momentum independent term uq has been explicitly written out in Eq. ([28]). In this Hamiltonian, the band 
energies e Pj of the original model are present in the new band energy and also the interaction term, and hence the shift 
Eq. (11) cannot be absorbed in the \i alone. Now the effective Hamiltonian explicitly satisfies the Shift theorem- (I) 
provided Uq also transforms under Eq. (11) as: 



^0 

Thus in using the auxiliary Fermi liquid we have the 



4 u t + uj. 



(29) 



• Shift theorem- (1. 1): An arbitrary shift of t and J, as in Eq. (12), can be absorbed by shifting uq, as in Eq. (29), 
and the chemical potential /x —> /x + u t . 




eff 



FIG. 2: The pseudo-potential W e ff in the momentum space representation. The zigzag line represents W e ff. Note that the 
momentum transfer in the argument of J is also expressible as J P1 - P4 - 

Since the standard notation for interaction reads ^(ab\V\a f b f ) f\fl fy f a ' for a conventional two body interac- 
tion, our notation corresponds to writing W e ff (pi,£>2 5P3,P4) = (piP2\W\p4Ps) . Fermi symmetry implies the in- 
variance W / e //(Pi 5 P2;P3 5 P4) = W e ff{p2,Pi\P4,P3), and Hermiticity implies the invariance W / e //(pi 5 P2;P3 5 P4) = 
W e f f(p3,p4]pi,p2). For this canonical theory, we calculate the commutator: 

[fiaii Heff] ~ ^ ^ ^ijfjcTi H> fieri H~ -^iai 



with 



3° '3 



[fi<Ti,V e ff]. 



9 Jij i^i^j) fjaifjvjfio 



(30) 



(31) 



Let us note that Eq. (31) differs from in Eq 



we impose the single occupancy constraint on the auxi 



20), through terms (in underbraces) that vanish identically if 



iary electrons. 



D. Equation of Motion for Q continued. 



We now return to the study of the equation of motion for Q in Eq. (24), given in terms of A ia of Eq. (20). Since 
this object represents the crucial Heisenberg equation of motion for the destruction operator, we look for terms that 
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can be added to it to make it identical to Eq. (31). Comparing Eg . (20) and Eq. (31) we see that these differ by 
terms (the second and third terms of the square bracket in Eq. (31)) that are automatically vanishing for the Xf 3 
operators on using their standard rules. Thus we can add such vanishing terms to Eq. (20) that remain exact and 
also importantly preserve the Hermitian nature of the auxiliary Fermionic theory in approximate schemes. We thus 
rewrite also an exact, but more useful result: 



13 a j 



1 

2" J 



1 



2 1 



Wja^X^X^, (32) 



so that Ai^ (Ti and A^ <Ti contain terms that are in one to one correspondence. We will use Eq. (32) in in place of 
Eq. (20) in Eq. (24) next. Let us recall also that Jij = J^- — S[ij] uq here. 

The notation simplifies if we use the matrix notation for the spin indices introduced in Ref. (10) and Ref. (|2j) e.g. 
Gcnaflh f] ~^ f]]ai<Tfi so that we may regard Q as a 2 x 2 matrix. In short, the space-time indices are displayed 
but the spin indices are hidden in the above matrix structure. We next define 7 through: 



[i or 



(33) 



where we denote the k conjugation of any matrix M by (M^) ai(72 = M^^aic^. This conjugation corresponds to 
time reversal in the spin space. Let 1 be the identity matrix in the 2x2 dimensional spin space. 
We employ a useful relation with an arbitrary operator Q that follows from Eq. (14): We write 



where we denote 



and 



la a a b [i,j] = {Wb) Ga b a a \jrr,i n } = ((<T a tr b Xf^Xf")), 



-±<r b L 



(34) 

(35) 
(36) 



3 5V- iaj 



OVA A 



andD[M] = D[i\. 



(37) 



In and we have are equal time objects with creation operators to the left of destruction operators. Let us 

note the rewriting of the last term in Eq. ( [24] ) : 

-E v S lff2 (^1) ((xr>(n)x^(n)xy°(T f ))) = +^V^7 w [i]-V 1 H)^ / b I /]- (38) 



With this preparation, using Eq. (32) and we rewrite Eq. (24) as 

(d T , -ti+\jo)g„ iaj [ij\ = -5[ij](5^, - 7W*D - vT'Vi) g^ f [ij\ - v°p{n) g^j] 

+V% a2 (7<™ (i) - D„ i<Tl (i)) g a2(7f (j, /) 

+ t[*,j]{ (i-7[i] + ^[i]-^7[i] + ^[i])-e[i,/] } + m]{(-\i\iA + \m,i])-g[iJ] 



(39) 



where the fixed variables are in normal letters and the (repeated) summed variables are in bold letters. This may be 
written compactly in matrix form as 

(d Ti - n)g[i,f\ = -s[ij](t - 7 [*]) - v, • g[ij] - Hj • g%f] + (7(0 - d[{\) • v 4J • g%f\ 

-x[i,j]-g\}j]-Y[i,j]-g\jj], (40) 
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where we used the definitions (with fixed j and summed k) 
X[i,j] -- 



*M (D[i\ + \d\i\) + S[i,j] \ (J[i,k] D[k] - t[ik] D[k,i] 



Y[i,j] = -t[i,j] (l-7[i]-\^[j}) + \j t-5[i,j]\{J[iM lM-t[iM 7 [k,i] 



(41) 



These exact equations Eq. (40) and Eq. (41) form the basis for the remaining discussion. The coefficients in X and 



Y differ slightly from the ones in Ref. (2), in view of the usage of the effective Hamiltonian idea in this paper. The 
extra terms arise from the form of Eq. (31 ), and actually vanish if we could treat either of these exactly. We will show 



that this formulation leads to approximations obeying the Shift Theorems (I-II) discussed earlier; note however that 
Eq. (40) and the forms of X,Y in Eq. (41) are manifestly invariant under these theorems. 



III. DECOMPOSITION OF Q INTO THE AUXILIARY FERMION GREENS FUNCTION g AND THE 

CAPARISON FACTOR fi 

As discussed in the Introduction, we next write the product ansatz for Q 

G[a,b] = g[a,r]./4r,6], (42) 

where g is the canonical auxiliary Greens function and \i is the caparison factor, or the adaptive spectral weight. Since 
Q satisfies antiperiodic boundary conditions under r a — > r a + j3 and 75 — » r& + /? separately, we must Fourier transform 
both factors g and \i with Fermionic frequencies uo n = (2n + l)7r/c#T. At this point \i and g are undetermined. Let 
us first note in matrix notation the equal time objects: 

7 [i] = a[r,i]^(g[i ) a]- M [a ) i])( fe ) = ( M [a ) i])( fe )-(g[i ) a])( fc ) ) 
7M = S\j~,i] -»• (gb>] -Ma,i]) (fe) = (M[a,i]) (fc) • (g[j,a]) (fc) . (43) 
We define a three point vertex functions 



or as an implicit matrix in the upper indices (but explicit in the lower ones): 



A<7 3 a 4 (p,g;r) = - ^ v f 3(74 {g 1 \p,q] }, W CT3(J4 [a,&;c] = ■ (45) 



In a similar vein, to obtain the four point vertex functions corresponding to the source V rs with a pair of points 
r, s, we define: 

A,(p,?;r-,s) = -^{g- 1 b,g]}, W.[a,6; C) rf] = ^M. (46) 

This set of vertices A and U replace the single vertex V of a canonical many body system, and we will also find 
equations determining these below. Clearly in any exact treatment, the four point vertex contains the three point 
vertex by collapsing the points: 

A££(p,g;r) = A££(p,g;r,* r), (47) 

and similarly for U. However in any approximation scheme, this identity would follow only if the single occupancy 
constraint at a given site i namely: (( X^ 1(T2 X® as . . .)) = is satisfied exactly, for all choices of the spin indices. Since 
typical approximations relax this constraint, if only slightly, it is therefore useful to keep both the sets of vertices in 
the theory as separate entities. Another attractive possibility is to require the identity Eq. (47), by making a different 
set of (controlled) approximations, and is also discussed below. Fig. Q illustrates the conventions used for the four 
point vertex, the three point vertex is obtained by the indicated contraction. 
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We now use a notation where * is used as a place holder, as illustrated in component form by 



'CF a CF h • • 



5V? a,ab ' 



with icr a (T h — &a&bi and an implicit spin flip in the indices of the attached derivative operator 5/5V? a,(7b . 
We would like to rewrite Eq. (40) in terms of the vertex functions. We need to express 



■ g% /] = -t[i,j] (£>[*] + \D[j}) ■ g[j, /] + l - k] D[k] - t[*k] D[k, i}) ■ g[i, /] 



in terms of the vertex functions. 



Differentiating Eq. (42) we find 

(5 



7 G[a,b] = g[a,c] • A CT5f7 /(c,d;r) • <?[d, 6] + g[a, c] • W^/ [c, 6; r] 



Consulting Eq. |37j ) for the definition of ^o- l0 - 2 [*] = (ai^ ^v^ 2 ? wnere — (J i (J 2^ we rewrite this as 

D[r]-g[a,b]=C -s[cl,c] - A*(c,d;r) • Q[d, b] + £* ■ g[a, c] • W*[c,6;r], 

where the spin flip in the derivatives is implied as stressed above. 

Combining Eq. (49) and Eq. (51) we define the useful linear operator 



S 1 S \ 1 SI. s 

W* + 2 5V*) + 2 t[iM C ' gM ■ SVl t ~ 2 J[i,k] C ' gM ' 5V*' 



Hence we may write Eq. (49) compactly as 

x[i,j].g%f] = *[*,b]-a[b, /] + *[*,/] 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



where the two central objects of this theory arise from the action of a common operator Eq. (52) on two seed objects 
g _1 and fx as follows: 



<E>[i, m] = L[i, c] • g 1 [c, m] 

-t[i,5\ r • g[j,c] • ( A*[c,ra;i] + ^A*[c,ra;j] 



^[i,k] <f • g[i,c] • A*[c,ra;k,i] + ^J[i,k] <f • g[i, c] • A*[c, m; k], 

(54) 



and 



^[i, m] = — L[i, c] • /i[c, m] 

— j] • g[j,c] • ^W»[c,m;i] + ^W»[c,m;j]^ - ^t[i,k] £* ■ g[i, c] • W,[c, m; k,i] + ^J[i,k] £* • g[i,c] • W»[c,m;k]. 

(55) 



Writing Eq. (41) as 



Y[i,j] = -*M+YiM 

= tM (7W + ^7b1)-^,i] i(J[i,k] 7 [k]-t[i,k] 7 [k,i]). 



(56) 



We also need to process the object: 



( 7 (») - d[*]) • Vij • /] = 7(0 • • g% /] - r • Vi j • ^ <?[j, /] 



7(0 • Vij • /] - r • Vij • g[j, c] • A,[c, r; i] • a[r, /]-£*• V 4J • g[j, c] • W.[c, /; i] 

(57) 
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IV. ASSEMBLING THE EQUATIONS 

Let us rewrite the three relevant equations symbolically: 



1. Eq. (40) for Q\ 



(d Tz - fi)g = -s(t - 7 ) - v l • g - Vi j • Q + (7 - A) • v. j • Q - x • a - y • a, 



2. Eq. (57) for the two site source Vij\ 



(58) 



( 7 - A) • Vij • 5 = 7(0 • ^ j • Q\j, f] - r • Vij • g • A, • - r • Vij • g • w. 



3. Eq. (53) the product rule: 



Combining these we rewrite Eq. ( [58] ) symbolically as 

(a ri -/x + y + V, + (l-7)-V,j+r - Vij - g-A* + ^)^ = -^(l-7)-^-r-V,j-g-^ 

Defining 

go 1 !*,/] = {[(M-Sr 4 -|Jb)i-ViMi ) /] + t[i ) /]-v i , / } ) 



the exact EOM Eq. (40) can be written in matrix form: 



(59) 



(60) 



(61) 



(62) 



(63) 



Clearly this becomes the exact equation Eq. (63) at A = 1, and reduces to the Fermi gas Greens function Eq. 



{go 1 [i.j]+7i-V i j-r-V i , a - g(a,b)-A,(b ) j;*)- Y 1 [i, j] - j]} • g[j, f] • M [f, /] = 
/] (1 - 7 W) + /] + T • Vi, a • g(a, b) • W,(b, /; *). 

At this point, a convenient parameter A (finally set A — >• 1) is now inserted into this equation as follows: 

{go 1 [i.j] + A7i-Vij-A£*-Vi,a- g(a ) b)-A,(b ) j;*)-A j] - A $[*, j]} • g[j, f] • /x[f, /] = 
/] (1 - A 7 [t]) + A /] + A T • V l>a • g(a, b) • W,(b, /; i) . 

s v ' 

(64) 
(|62| at 

(65) 
(66) 

We can usefully invert Eq. ( [65] ) and write 

g- 1 [i, m] = {g^ 1 [i, m] + A 7i • V», m - Af ■ V <>a ■ g(a, b) • A*(b, m; i) - A Y 1 [i, m] - A m]}. (67) 

We see that g satisfies a canonical equation, with a delta function of weight unity on the right, and /i soaks up the 
remaining factors on the right hand side of Eq. (64). This is decomposition is not unique, one has the obvious freedom 



A = 0. We may now split Eq. ([63]) exactly into a pair of equations that are fundamental to the theory: 

{go 1 M+A7i-V<j-Ae*-V< 1 a- g(a,b)-A*(b,j;z)-A Yx[i, j] - A j]} • g[j, /] = 5[i, f] 

/] = 5[i, f] (1 - A 7 [i]) + A /] + A T • V,, a • g(a, b) • W*(b, /; z). 



of respectively post-multiplying g and pre-multiplying \i by a common function and its inverse. However, requiring g 
to be canonical fixes the function to be unity. The motivation of introducing A in the above equations, is to establish 
adiabatic, or more properly, parametric continuity with the Fermi gas 14 . At this stage some remarks are necessary 



• At A = 1 Eq. (65) and Eq. (66) becomes the exact equations for the EC phase, while it has the virtue that as 



A = it gives a canonical equation for g, with j] = 15[i,j]. Procedurally, we can calculate objects to a given 
order in A iteratively, and set A = 1 at the end of the calculation. We thus establish and maintain continuity 
with the Fermi gas in the equations of motion. 
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The process of introducing A into the EOM is not unique. For example the terms of Eq. (64) in the underbraces 



cancel at i = j from the vanishing of Eq. (38). However this cancellation is exact only at A = 1, so we will find 
below that an expansion in A has the annoying feature of a slight violation of the contraction of indices result 
Eq. (47). We will show below that this is inconsequential to the orders in A considered here. With hindsight, a 
better strategy would be to impose the constraint Eq. p7| ) to the order of the calculation. This can be achieved 
if we multiply the terms in underbraces by a sufficiently high power of A r , say with r > r$ , and thereby avoid 
dealing with this problem at low orders r < tq. Below we will analyze the minimal choice r = 1, record the 
issues that crop up and make suitable approximations later. The impatient may simply ignore the terms with 
underbraces. 



Another type of freedom is available at this stage: if necessary, we could add an arbitrary term that varies 
smoothly with A and vanishes at both end points e.g. oc A(l — A), to either side of Eq. (65) and Eq. (66). It will 
turn out that the first order term [ g _1 ]i calculated below, does need a simple term of this type to fulfill the 
Fermi surface sum rule. In general, however, the natural and minimal choice made in Eq. (64), without such a 
term, seems adequate for higher terms. 



• We note that the Shift theorems (I-II) are preserved by X, Y above in Eq. (41 ), and this invariance survives the 
introduction of A in Eq. (64). As a result the various objects <£, g _1 , \i satisfy these theorems individually. 



This property leads to a powerful consistency check on the approximations to each order in A. 



Note that a A expansion of implies that the high frequency fall off of the Q ~ gj, now occurs with a coefficient 
Co = 1 — A 7 that is different from 1 — ^ at finite orders of A. While it is tempting to freeze this coefficient at the 
exact value, it would be inconsistent since we take its derivatives to find ^ etc. The departure of this coefficient 
from the exact value becomes increasingly significant near n ~ 1, and provides a criterion for the validity of a 
given order of approximation. 



V. EXPLICIT EQUATIONS AND THE ZERO SOURCE LIMIT IN FOURIER SPACE 

When we turn off the sources, the various matrix function Q, g, /i become spin diagonal. We will also take Fourier 
transforms (only) in this limit, since translation invariance in space and time is regained when the sources vanish. 
We next express and \I/ explicitly in terms of the vertex functions. We need to take the Fourier transform of 



Eq. (54) and Eq. (55). In the ECFL theory, a rotationally invariant liquid phase is obtained by turning off the sources. 
We can use the standard spin rotational symmetry analysis illustrated here with A as in Ref. (fTUj) . We define the 
three non vanishing matrix elements as A^) = A££, A^ 2 ^ = A££ and A^ 3 ^ = A££. We also record the Nozieres identity 
for the two expressions of a particle hole singlet: A^) — A^ 2 ^ = A^ 3 \ which provides an important check on the theory. 
We further use a notation for the frequently occurring antisymmetric combination A^ a ^ = A^ 2 ^ — A^ 3 ^. Armed with 
these, we next drop the matrix structure by utilizing an identity arising with a fixed a (such as in the expression for 
$> a(T above): 



<<r| £* • g • A* \a) = Y, ™* Sa a a b Kf a 

— (gcrcr — Sad- ^ao) 

= g(A^ -A^)= gA^. (68) 

Note that we dropped the spin index on g due to the isotropy of the state. 

We use the FT convention for the two, three and four site objects illustrated with the examples: 

g[ a ,b] = J2 eik{a ~ b) $( k ) 

k 

Pl,P2 

A%£l[a,b;c,d\ = £ e i ^ a -^ b -^ c+ ^A^l(p 1 ,p 2 -,p 3 ,p i ). (69) 

Pl+P4=P2+P3 
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The identity Eq. (47) in momentum space implies: 

A(pi,p 2 ) = Yl A (^i^2;P3,P4), 

P3,P4 



(70) 






^4 



0-30-4 L 

and 

Kl2\PhP2]P3]P4] 



FIG. 3: The vertex The four site vertex can be visualized from its definition for canonical theory: —{{f a a 1 fl a2 fl a 3 fd <t 4 )) + 

{{fa cr-tfl cr 2 )){{fc (T 3 fd cr 4 )) = §v cr 3 a 4 gcr lC r 2 (pb) — g CT1CT / [a, a] A^ j a> b; c? g , a2 [b, 6]. Therefore we may visualize that apart 

cd 

from the external legs, &%\% 2 A [a, 6; c, d] ~ (/ a ^ifl a2 ft a 3 fd <x 4 )- Note that in this convention, the labels differ by a cyclic 
permutation from those in Fig. 



At zero source we get the exact relations between self energies and vertices by Fourier transforming Eq. (54) and 
Eq. ([55 ) 



*( fc ) = ^\£p + \ e k + \jk-p\ g[p] A (a) (P,^) + ^^g+p-fc g[p] A {a) (p,k;q+p-k,q) 

P \ / pgr 

= X) + ^ £/c + \ £ Q+P-k + ^-p) SW A ia \p,k;q+p- k,q) 

pq ^ ' 

= ( £ p + + ^ e «+p-fc + ^ J fc-p ) gW U {a \p,k]q+p-k,q) (71) 



A convergence factor e zu;p0+ arises from the time ordering and is implied wherever necessary and the last line in both 
equations is valid provided the identity Eq. (FFob is satisfied. Here A^ a ) = A^ 2 ) - A^ 3 ) and = - . 



With /c = (fc, zcjfe) and co> n = 7r(2n + l)/c#T, the Greens functions at a fixed A read: 

g(k) = g (fc) x fi(k) 

g _1 (fc) = iw n + /x-£ fe - ij - A Fi(fc) - Ai>(fc) 
fi(k) = I-A7 + A 



(72) 
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The sum rule for the number of physical particles and the auxiliary Fermions is given by 

n 
2 



p 



n 
2 



(73) 
(74) 



While the sum rule Eq. (87) clearly counts the number of physical electrons, the origin the sum rule Eq. (74) for g 



requires some discussion taken from Ref. ([2]). We recall that it is meant to enforce the Luttinger Ward theorem of 
a conserved Fermi volume for the auxiliary Fermions. By so doing and through the composition Q = g x /i, it also 
preserves it for the physical Fermions. While \x provides us with one obvious Lagrange multiplier to enforce one 
of the sum rules, the more subtle parameter u Q: introduced in Eq. p5|), is required to enforce the second sum rule 



Eq. (88). Explicit expressions for 7 , Yi,<I>,\£ can be calculated order by in A as demonstrated below. 



VI. SUMMARIZING 



Before proceeding to the iterative scheme, we collect all the relevant equations for convenience. 



Table 1: Summary of all Equation 



Object 


Defining Equation 


Eq. 
No. 


g _1 [z,ra] 


{go 1 [« J 77i] + A7i-V i , m -Ar-Vi,a- g(a,b) • A*(b,ra;z) - A Y\[i, m] — X 3>[z, m]} 


Eq. ( 


65 
66 
56 
43 
43 
54 




H[i,m] 


5[i,m] (1 - A 7 [z]) + A ^,m]+Af • V i>a ■ g(a, b) • W*(b, m; z) 


Eq. ( 


Yi[z,ra] 


t[i, m] ( 7 [z] + | 7 [m]) - m] ± ( J[z, k] 7 [k] - t[i, k] 7 [k, z]) 


Eq. ( 




/i (/c) [a, z] • g (fe) [z,a] 


Eq. ( 


7 [z,ra] 


/i^[a,z] • g (fe) [m,a] 


Eq. ( 


<3>[z, m] 


r • g[j,c] • (A*[c,ra;z] + ±A*[c,ra;j]) - k] f • g[z,c] • A*[c,ra;k,z] + 
|j[z,k] <f -g[z,c] • A*[c,ra;k] 


Eq. ( 


\I/[z, m] 


-t[hi] r • g[j,c] • (W*[c,ra;z] + §W*[c,ra;j]) - §t[i,k] f* • g[z,c] ■ W*[c,ra;k,i] + 
|j[z,k] C -g[z,c] -W*[c,ra;k] 


Eq. ( 


55 


) 



The various vertex functions are found from: 

Table 2: Vertex functions 



Vertex 


Defining Equation 


Eq. 
No. 


A°*J>[i,m;j] 


(sV CTcCTd ) £°~aO- b m ] 


Eq. (44) 


A%' e %[i,m;j,k] 


( sv l~ d ) g^M 


Eq. (46) 


K:^[i,m;j] 


( 5 y CT c^ d ) f^CTaCrblh 171 ] 


Eq. (44) 


UZ£i[i,m;j,k] 


3,k 


Eq. (46) 



It is worthwhile providing one non trivial example of the matrix notation. In component form note that 3>[z,ra] 
can be written out as: 

$^ a Jz,m] = -t[z,j] Giox g CTlCT2 [j,c] ^A^[c,m;z] + ^A^[c,m;j]^ 

- ^t[h k] (JiCJ! g ai(T2 [z, c] A%?™ [c, m; k, z] + ^ J[z, k] g ai(T2 [z, c] A^ m [c, m; k], 



VII. A EXPANSION AND THE ITERATIVE SCHEME 



Taking functional derivatives w.r.t. V, we generate a self energy - vertex hierarchy of Fermionic theory, paralleling 
the standard (i.e. canonical) theory, but with greater complexity due to the two kinds of vertex functions and self 
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energies. We describe the A expansion and the iterative process next. The iterations are analogous to the skeleton 
diagram expansion in standard many body theory, where A plays the role of the interaction constant. Various objects 
are expanded in terms of A and g, while g itself is left intact. Potentially confusing is the treatment of g _1 , which 
is expanded in A and g, ignoring its obvious relationship as the inverse of g. This becomes understandable when we 
recall that g _1 is, apart from g^" 1 , the Dyson self energy of the system, and is to be regarded as a functional of g, as 
in the Luttinger Ward functional Ref. (17). One example of this expansion may be useful. Consider 7(2, ra], we will 
expand it as: 



7 (/c) [i,m] =g[m,a] -/i[a,i] =g[m,a] • ([/x[a,i]] + A (/i[a, i}} 1 + A 2 [/x[a, i}} 2 + 0(A 3 )) 



(75) 



keeping g intact, i.e. unexpanded in A. A similar expansion is carried out also for leading to a correction of the 
high frequency fall of coefficient Co as noted above. 

Iterative process: We now describe the various steps of the iteration process. First note that all variables (except 
g) are expanded as 



A = [A]o + \ [4]i + A 2 [A} 2 - 



\? [A] T 



(76) 



The iteration scheme can be summarized in the two following tables. Table 3 lists the seed objects needed at any 
order and gives the derived objects. 



Table 3: Level p calculations 



Seed object 


Derived objects 


[ fi[i,m] } p 


py[z], 7[i,m], Yi[i,ra], W[a,6;c], U[a,b\c,d]] p 


[ W[a,6;c], U\a,b\c,d\ } p 




[g~ 1 [hm]] p 


[ A[a,6;c], A[a,6;c,d] ] p 


[ A[a,6;c], A[a,6;c,d] } p 


Mhm]] p 



Table 4 lists the higher order objects and the needed lower level objects for stepping up. 

Table 4: Level step-up calculations 



Level (p + 1) object 


Required level p objects 


[ fi[i,m] ] p+ i 


[7[i], ^[i,m], W[a,ft;c], W[a,6;c,d]] p 


[g _1 M] p+1 


[Yi[z,m], $[i,ra], A[a,ft;c] ] p 



I. Initialization at p=0: The iterations require the following starting relations. 

go^i.m] = {[(p-d Ti -^J )l-Vi]5[i,m]+t[i,m]-Vi, m } 
Hi,f}} = 1 *[*,/] 



(77) 



• II. Computation of derived objects at level p from Table 1 : 

The set of equations requiring \p[i,m]] 



h[i\} P 
h[h m ]] p 



•g (fc) M 



■g (fc) [m,a] 



'£$[i,m;j]] p = (tt^) K^Mlp 



K:Z»[i,m;j,k] 



( Z\Wc<Jd ) [Mo 

V j,k 



Yi[i,m] ] 



t\i, ml 



7[iJ + -7[mJ 



■*[i,m] ^[J[i,k] 7 [k]-i[i,k] 7 [k,il 



(78) 
(79) 
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[ V[i,m] } p = -t[i,i] T • g[j,c] • ^Z4[c,ra;i] + ^W*[c,m;j] 



^[i,k] r-g[^c]- (W*[c,ra;k,z] ) p + |j[i,k] f'gfrc]- (W*[c,ra;k]) ? 



T/ie se£ of equations requiring [g ^ra]] 



[ A^[i,m;j] ] p = -(^^j) [g<^M ] 



[A^[i,m;j,fc] 



[ $[i,m] ] = -t[i,j] C ■ g[j,c] • A,[c,m;i] + -A,[c,m;j] 



-i[i,k] £*-g[i,c]- (A»[c,m;k,i] ) p + rJ[i,k] £*-g[i,c]- (A»[c, m; k] ) p , 



(80) 



(81) 



• III. Level p to Level (p+ 1): step up equations: 

H*."»]]p+i = -<*[*, m][7[»]]p+ [*[*,"»]]„ + [r-Vi, a -g(a,b)-W,(b,m;i)] p , 
[g-^m]]^ = [ 7i • Vi, m - T • Vi, a • g(a,b)-A,(b,m;i)] p -[yi[i,m]+*[i,m]] I , 

• IV. If required level is reached exit, else return to Step II. 



(82) 



Having formulated the iterative process, we next apply this to obtain the second order Greens functions. The 
calculations are detailed in the Appendix [Bj and we directly present the first and second ord er re sults here. Disp laying 
the so far hidden uq coefficient, we write the complete set of equations to 0(A 2 ) from Eq. (180) and Eq. (178). 



n[k] 

[$(*)]l 



2 

1 - A J + A 2 ^- - A 2 ( £ p + £k+ q -p + s k +£ q + Jk- P - uq) g[p] gfa] gfe + k - p] + 0(A 3 ) (83) 



iuo n + ia — 1 — A n + A 



3n- 



]^ + ^ \ J *-« S[Q] ~ A 2 [$ (*)li + 0(A 3 ) 



g[<?] g[p] g[fc + 9-p] 



(84) 



X (gjfe + £ p + £ g + Sk+q-p + Jfc-p — ^0 ) {?k + £p + £q + ^k+q-p + ^ (^-P + ^P-<?) — ^0 }• (85) 



Eq. (84) clearly satisfies the S hift theorem (I); it also satisfies Shift theorem (II), provided we restore the omitted 
constant terms from Eq. (177). In using this expansion, one must first set A — >• 1. The self energy from a Feynman 



diagram theory to second order from H e fj m Eq. (25) matches the above expression for g _1 . The required diagrams 
are shown in Fig. Q up to second order where the zigzag line W e ff is defined in Fig. (J2|. Apart from a single term 
(the expansion of Y\ in A), the expansion of the auxiliary Fermi liquid is largely "autonomous", i.e. proceeds without 
requiring the knowledge of /i, and is represented in Feynman diagrammatic terms. The caparison term \i has no 
obvious interpretation in terms of H e ffi but is easy to compute along lines similar to the ones shown here, and the 
full theory ties up the two factors in Q in a specific way as described here. 

A consistent first order, i.e. 0(A) theory for g _1 and \i can be found after dropping all 0(A 2 ) terms. As it stands, 
we would get /i = 1 — X^ to this order, and this would violate the Fermi surface volume theorem ( Ref. (13)). We may, 
however, also set ji[k] to unity. Formally this is achieved by adding A(l — X)^ to ja[k] as discussed below Eq. (64), 
since this added term vanishes at both limits of A. This procedure is within the permissible adjustments of the 
continuity argument, and at second order cancels out so that the quoted second order result is unchanged. Further all 
vertices are unchanged since this is a static term. In this way the first order theory can also be arranged to satisfy the 
Luttinger Ward Fermi volume theorem. This theory has a ba nd dispersion (1 — n)ek that that shrinks in width by a 
factor (1 — n) as in the Gutzwiller-Brinkman-Rice theor^^, with an enhanced effective mass m/m* = (1 — n). The 
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cr a rj & 



FIG. 4: The self energy graphs to second order from W e ff and the effective Hamiltonian H e ff. These determine the $ self 
energy. 



second order result here corrects the above result in a crucially distinct frequency dependent fashion, as discussed in 
this paper. 

The physical Greens function is obtained from: 

g[k] = s [k]»[k]. (86) 

The number of the physical electrons is fixed by the first sum rule: 

\ = £e[fc]e^ 0+ , (87) 

k 

while the auxiliary Fermion satisfy an identical sum rule: 

\ = EsM eiaj " 0+ ( 88 ) 

fc 

We can determine the two independent real parameters fj, and uq in order to satisfy both these equations simulta- 
neously, and thus the role of uo as a Lagrange multiplier, similar to that of \i is now evident. It is also clear that 
the shifts of t or J can be absorbed in the two Lagrange multipliers fi and uq. It is worth noting that the Simplified 
ECFL model used in Ref. (|2j) and Ref. (|8j) can be obtained from Eq. (83) and Eq. (85) by throwing out the band 
energies and exchange energies in the coefficients of g(#)g(p)g(fc + q — p) while retaining uo, so that the Lagrange 
multiplier of that approximation Ao, is essentially the same as u$. 

The role of the two sum rules in fixing the number of Fermions and also the Luttinger Ward Fermi surface is already 
discussed in Ref. (|2]). We can add to that discussion with the help of the explicit functional forms found above. It 
should be noted from Eq. (83) and Eq. ([85]) that the functional derivatives 



are symmetric functions under k p. This symmetry therefore guarantees the existence of two Luttinger Ward type 
functionals of the auxiliary Greens function g, 

Oi>[g] = ^ W(k,q\r,p) [W(k,q-,r,p) + W(k,q;p,r)} g[fc] g[p] g[g] g[r] 

k,p,q,r 

= J E W(k,q;r,p)s[k] S \p]s[q]s[r], (90) 

k,p,q,r 

such that the two self energies can be found from these functionals: 



(91) 
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The form of these two functional follows to this order from Eq. (85), and it is natural to conjecture that such 
functionals exist to all orders in A. The existence of the fl$> functional guarantees a (FS) volume conserving Luttinger 
Ward Fermi surface for the g electrons, and the smooth behavior of ^f(k) near this surface guarantees likewise for the 
physical electrons. 

VIII. WARD IDENTITIES 

This theory admits Ward identities involving the vertices A and U that guarantee current conservation in a similar 
fashion as Ref. (fT0|). This is displayed with the help of sources, the charge potential u[m] = ^ a V^ 7 and an added 
source v[m] coupling to the kinetic energy as 

t[i,j]^t[i,j](l + v\j]-v[i\), (92) 

so that v[j] — v[i] acts as a discrete version of the Peierls phase factor of electromagnetic coupling in tight binding 
systems. We define 

5 5 

D m =d Trn -—-—- (93) 
ou[m\ ov[m\ 

so that the Ward identity expressing the conservation of current, from Ref. (10) reads 

An G[iJ] = (5[i,m]-S\fM) G[hf\- (94) 

This is a discrete (Takahashi type) version of the usual Ward identity appropriate to the lattice Fermi system at hand, 
and electromagnetic coupling only requires the long wavelength limit of this identity. We will define the (T) vertices 
(summing over a) 

5 

A {T) (i,j;m) = -j^g-*[i,j]/ U:V ^o 

U ir) (iJ;m) = ——-{i aa [ij}/ u , v ^ . (95) 
ov[m\ 



It is easy to see that the bare r vertices are given by differentiating g 1 in Eq. ( 62 ) as 



A (r) (i,j;m) =t[ij](5[i,m] -6\j,m]), X {r) [p u p 2 } = e Pl - e P2 , (96) 

while the singlet (i.e. density) vertices are already known from A^ s ^ = ^2 a(jf ^V'a" Note that the (T) type vertices 
are antisymmetric in i O j or p\ O p 2 . 



Taking Fourier transforms in Eq. (94) and writing Q = g x /i, we get the conservation law: 

{iU) pi -iu p2 ) (g(pi)A (s) (Pl^2)g(p2)/i(p 2 )+g(pi)^ (s) (Pl^2)) - (g(Pl)A (T) (Pl^2)g(p2)/i(p2)+g(pi)^ (T) (Pl^2)) = 

g(P2)/i(P2) - g(PiMPi)- (97) 

Canceling out g(pi)g(p2) we get the Ward identity: 

(i0J pi - i0J p2 ) ( y A {S \p 1 ,P2)KP2)+H {s) (pi,P2)E- 1 (P2)) ~ (A {T) (pl,P2)n(p2)+U {T) (pi,P2)g- 1 (p2)) = 

g" 1 (Pi)^(P2)-g- 1 (p 2 )M(Pi)- (98) 

With ito n —> z n , we rewrite this as 

Wg(Pl,P2) /i(p2) + g~ 1 (P2) W^(PUP2) = 0, (99) 

where we have defined the two Ward functions: 

Wfl(pi,P2) = (zi - z 2 ) aW( Pi , P2 ) - ^ r H Pl ,P2) + g- 1 ( P 2) - g-Hpi) 

W P (pi,p 2 ) = (zi - z 2 ) U (s) { Pl , P 2) - U {r \ Pl ,p2) + M(Pi) - M(P2) (100) 
Since p\ and p 2 are arbitrary, the two terms must vanish separately giving us the pair of Ward identities: 

W fl (pi,P2) = 0, (101) 
W M (pi,P2) = 0. (102) 
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IX. RANDOM PHASE APPROXIMATION 



Since the Greens functions are known to 0(A 2 ), we can take the derivatives of Eq. (137) and Eq. (138), to get 



vertices to this order. Here we calculate by taking the equations to 0(A) only, but assuming ^g = gAg rather than 
^g = gg? thereby obtaining the analog of the RPA. Since the spin susceptibility is also of considerable interest, 
we will calculate the required vertices in the that channel as well. Summarizing the results we write linear integral 
equations for the U vertices: 

U (r) \pi,P2] = -A^g( 9 )A( r )( (/ , (Z + p 2 -p 1 )g( (/ + p 2 -p 1 ) + 0(A 2 ) 
q 

U (s) \pi,P2] = -A ^ g(g)AW ( 9 , 9 + P2 - Pl )g( g + p 2 - Pi ) + 0(A 2 ) 

U (t) lPi,P2] = A^g( g )A(*)( g , 9 + p 2 -p 1 )g( 9 + p 2 -p 1 ) + 0(A 2 ), (103) 



and similarly for the A vertices: 

A (T) [pi,p 2 ] = (e Pl -e P2 ) (l- Aji) - A ^ g(g)A (r) (<?, q + p 2 - pi)g(q + p 2 - Pi) J r {q,Pi,p 2 ) + 0(A 2 ), 

q 

A (s) bi,P2] = 1- \^2e(q)A ( -'\q,q + p 2 -pi)e{Q+P2-Pi) ^(<Z,Pi,P 2 ) + 0(A 2 ), 

A (t) bi,P2] = i + xJ2s(q)^ {t) (Q,q + P2-Pi)s(q + P2-Pi) J r (q,Pi,P2) + o(\ 2 ). 



(104) 



where we use the shorthand J r (q,Pi,p 2 ) = {z Pl + £ P2 + £q + £ q+p 2 -pi ~ u o + \ (Jpi-p 2 + Jq-pi)}- The term in 
underbrace receives an 0(A) contr ibuti on from differentiating the explicit v dependence of the transformed t[i,m] — > 
t[i, m](l + v[m] —v[i]) term in Eq. ( 137). It is readily shown by examining the kernel of the integral equations that the 
solution for (pi,p 2 ) is antisymmetric under exchanging pi p 2: while A( s \pi,p 2 ) and (pi,p 2 ) are symmetric. 

These vertices are shown to be compatible with Ward identities to 0(A) if used with the first order versions of the 
Greens functions Eq. ([83]) and Eq. ([84]): 



A 



g(p) = iuj n + /x - (1 - An)^ + - ^ J/e-g g(g) + 0(A 2 ), and = 1, 



(105) 



by substituting in the expressions Eq. (101 ) and Eq. (102), and showing the self consistency of this result. The details 



of this verification parallel the standard proof in QED and are omitted here. Note that \i must be chosen to be unity 
rather than 1 — as discussed in the second para below Eq. (84), although this choice is irrelevant to the verification 
of the Ward identity. 



X. TWO PARTICLE RESPONSE 



We are interested in the pair correlations of the density n a 
of the Greens function 



X° a and the spin density S* 



/ T,f _ X^ 1<72 , where r z is the usual Pauli matrix. These can be obtained from taking the functional derivatives 



(106) 



and can be conveniently found from taking a limit of the three site object Y£i£ 2 (p, q; r) = 8V ^ A Q ai a 2 [P? With the 
singlet and triplet objects denoted with a superscript a = s, £, we note the following relationships with the standard 
charge and spin susceptibilities of interest: 



((n a (r a )n b (r b ))) = n 2 -2T^(a,b) 



(107) 
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Owing to the Bosonic nature of the densities, we have the symmetry T^ a \b, a) = Y( a \a,b) from which the Fourier 
transform at Q = (Q,iQ q ) satisfies the relation: 



T< a >(Q) = T< a >(-Q). 

This symmetry can be used as another test of the consistency of any approximation. 

The Greens function in Eq. ( 106 ) can be decomposed in to g and \i as before and we find 



(108) 



{g<r lCTo (a,a)/v a<T2 (a, &)} , 



8Vr 3<74 

{g(a,b)A ff3CT4 (b,s;r)g(s,aV(a,q)} + {g(a,b) £4 3<T4 (b,&; r)} 



(109) 



where the vertex and T carry upper spin indices that are part of the matrix product. Turning off the sources, we find 
the expressions for singlet and triplet response 



T^\a,b;r) = g(«,b)A( a )(b,s)g(s,aHa,6) + g( a ,b)W tt (b,6;r), 

r (a Hpi,P2) = g(Pi)A (a) (pi,p 2 )g(p 2 )MP2) + g(Pi) u {a) (pi,P2), 



(110) 



where a = s, t. The definitional distinction between left and right derivatives leads to the asymmetry in the above 
equations making it necessary to test the consistency Eq. (108) term by term. 
Using the zero source limit notation from Ref. ([T0j) : 



= QZ; 
r>(°) - r>( 2 ) _ 



0(3). n {s) _ n (l) , n (2). 



)(«) _ 



(2) = Q(3) 



(111) 



The charge a = s and spin a = t susceptibilities at finite Q = (Q, are given by setting p 2 —> p and p\ p + Q 
and summing over p. 

T (a) (0) = E T(a) (P>P + 3) = E (s(p) A^(p,p + Q)g(p + Q)Mp + Q) + g(p)w (a) fep + Q)), (H2) 



These are exact expression for the susceptibilities, but as usual require a knowledge of t he ve rtices and Greens functions 
to give practical results. We can now use the RPA vertices calculated in Section. (IX) to give the corresponding 
expressions. 

We denote the susceptibility of the auxiliary Fermions as 



xi a) (Q) = -J2 S(q) A (a) (q, q + Q)g(« + Q), 



(113) 



and within RPA we note that is independent of p, and from Eq. (103) we denote that the U vertices are functions 
of the momentum difference only: 



U (a) \pi,P2} =\tia X { K \P2-Pl) 



(114) 



where £ a is 1 for a — singlet and — 1 for a = triplet. Therefore we can sum over the p dependence of the second term 
and rewrite Eq. (112) as 



(r^(Q)) RpA = -c aX ( «\Q) 



(115) 



where C a = (/i — £ a A §). It seems more appropriate to reset /i = (1 — A|) from unity at this level, in order to 
recover the expected high frequency behavior in the charge as well as spin channel, so that C S i ng i et = (1 — An) — > 1 — n 
and C triplet = 1- The vertices A are to be computed from Eq. (104) and form a consistent set of equations for two 
particle response in the sense of the usual RPA. 

The integral equations must be solved numerically. However in order to display some flavor of the results, we 
pursue this to the lowest order in A by iteration, where explicit results can be obtained. Let us define a few frequently 
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occurring generalized polarizability functions for convenience. We will now reinstate Jfe — > Jk — Wo 

Xo(Q) = ~^2g(q) g(q + Q) 

q 

xi(Q) = ~Y1 s(q)s(q + Q) {e q + e q+Q } 

q 

X2(Q) = l^2g(r)g(r + Q)g(p)g(p + Q) J P -r 
F(p + Q,p) = J2g(r)g(r + Q) |e p + e p+Q + e r + e r+Q + ^(J Q + J p _ r )| . 



(116) 



Here Xo(Q) is the standard Lindhard function and is positive in the static limit as Q — >• 0, while the other functions 
are generalizations thereof. 
The answers are 

T<'>(Q) = -(1 - An)xo(Q) - A[2 X o(Q) Xi(Q) ~ («o - ^o) X§(Q) + X 2 (Q)] 

T W (Q) = -Xo(Q) + A[2 X0 (Q) X i(Q) + (^Jq - «o ) Xo(Q) + X2(Q)]- (117) 
It is clear that the role of uq enhances the spin susceptibility while decreasing the charge susceptibility. To this order 



we see that the parity test Eq. (108) is satisfied to this order by using the symmetries of the objects in Eq. (116). 

Since the Greens function remain infinitely sharp within the RPA, its usefulness is limited- especially in view of the 
large frequency dependent corrections with characteristic asymmetry seen in second order results Ref. (|2j), Ref. J5J) 
and Ref. (9). A second order version of RPA seems most desirable, although even without vertex corrections to 
second order, the single particle spectral results are very interesting already. It also seems interesting to also study 
phenomenologically, the analog of the "bubble" diagram for purposes of extracting the optical conductivity; a scheme 



that reflects the width of the physical Greens function and satisfies the parity requirement Eq. ( 108 ) is given by: 

P{Q)\ ph en = ~Y^/2 £ G{q)G{q + 0) ' (118) 



although this expression is not the result a systematic expansion of Eq. (112) . 



XI. DISCUSSION AND CONCLUSIONS 

We have described above a controlled technique of dealing with the t-J model. This extremely correlated Fermi 
liquid theory is a strong coupling approach, specifically designed to deal with a hard many body problem. The 
considerations begin with the strong coupling limit of the Hubbard model, leading to the t-J model with a hard 
constraint of eliminated double occupancy. The Schwinger method gives us a crucial initial platform to deal with 
this problem. The ensuing exact functional differential equations are made tractable by the introduction of the exact 
product ansatz: Q = g x /i, with g a canonical Greens function of auxiliary electrons and /i the caparison factor. The 
latter, in turn, is understood as an adaptive spectral weight balancing the requirements at the high and low frequency 
ends of the spectrum. Both objects are expanded in powers of a parameter A, that plays the role of fractional double 
occupancy. Thus A = 1 corresponds to complete elimination of double occupancy whereas A < 1 has some residual 
double occupancy. We thus replace the hard constraint: of complete elimination of double occupancy by a softer one 
or partial removal. In order to provide a natural description of the canonical electrons, we introduce the effective 
Hamiltonian H e ff, depending parametrically on A. In order to obey the Shift theorems (I-II), we find it obligatory 
to (re) introduce a Hubbard type uq parameter in this model. It also plays the role of a second chemical potential 
as explained above. The set of steps followed, in our starting as well as ending up with a Hubbard type interaction 
has a slightly circular feel to it. This recipe is perhaps best understood as a renormalization group type procedure, 
where the constraint of single occupancy is enforced incrementally and the density of doubly occupied sites is thinned 
out smoothly. The infinite starting value of U in the t-J model is pushed downward to Uq, typically a fraction 
of the bandwidth from our numerical studies, albeit in a more general model H e ff, and is therefore amenable to 
a perturbative expansion. The form of the H e ff and the important role of the shift symmetries in validating the 
approximations is noteworthy. The hopping tij is elevated to an interaction constant of the model, this unfamiliar 
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step is kept under check by requiring the two important shift invariances. The Schwinger equation Eq. (40) for Q, 



being an exact statement of the problem, provides us with a rigorous backdrop to the entire procedure. Further our 
procedure has the advantage of being systematically improvable through the iterative scheme developed here. 

We can explore superconductivity at a qualitative level, by studying the pairing instabilities of the auxiliary Fermions 
given by H e ff via its BCS gap function A(fc). In this first approximation, the physical electron order parameter 
(Xfc°X^° k ) is proportional to that of the auxiliary electrons (/|(fc)/|(— fc)), together with the single occupancy con- 
straint of vanishing upon summing over the wavevector k. Within a generalized Hartree Fock theory, retaining the 
self energy correction to first order (as in Eq. (105)) as well as the pairing field average, we obtain an equation for 
the gap function A(fc): 



A(fc) 



1 



^0 



Jk-p \ A(p) 



tanh/3S(p)/2 
W) 



(119) 



fj,. In the computation below, we will neglect the 



where E{p) = yj A 2 (p) + and £ p = e p (l - n) - \ J q - P i 
numerically small J term in the single particle energy. Other than uq and the two single particle energies in Eq. (119) 



required for satisfying the Shift theorems (I-II), this is the same equation as the one found within the resonating 
valence bond theory in Ref. (|3J), Ref. ([15]) and Ref. (19). The tra nsitio n temperature for a d-wave state with a gap 
function A(fc) = A^ [cos(fc £C ) — cos(k y )] is obtained by solving Eq. (119) for the case of the nearest neighbour square 
lattice t-J model, with parameters indicated in the caption. It is straightforward to see that the T c equation has a 

250 



200 
150 
100 
50 




0.0 0.1 




FIG. 5: The transition temperature in Kelvin, from solving Eq. (119) assuming t — 3000K and J — 900K. The solid line 



indicates the likely regime of validity of the O(A) theory. Its dotted extension to lower hole density is speculative and is most 
likely to change with higher order corrections reflecting the nearby Mott insulating state. The dotted red line indicates the 
maximum T c obtainable from this scheme, and is seen to depend solely upon the magnitude of J. 



maximum scale of order J/(4&b) as already noted in Ref. ([15]) and Ref. ([T9]) . This value is attained in this solution 
at a higher particle density, or equivalent ly, a lower hole density, than is warranted by the first approximation. The 
solid line represents a plausible regime of validity of this scheme. 

The extended s-wave order is usually described by a gap function A(fc) = A S5 o + A Sj i [cos(k x ) + cos(k y )]. The 
constant term A S) o leads to a finite probability of double occupancy, since it survives a wavevector sum. After it is 
dropped as p er th e above discussion, the assumed (purely extended) s-wave order is supported by the J term in the 
kernel of Eq. (119) , but not by the uq dependent and single particle energy terms. The latter thus do not play a role 



in determining T c for either d-wave or s-wave orders despite their large magnitude relative to J. 

A detailed calculation of the gap equation is planned for the pairing of physical particles, parallel to the 0(A 2 ) 
theory of the normal state. The finite lifetime effects are then expected to become relevant. Such an improvement 
of the pairing scheme should yield a greater understanding of the balance between the different orders and a greater 
range of validity in density than the schematic theory treated here. 



XII. ACKNOWLEDGEMENTS 



I have benefitted from discussions with T. Banks, G.-H. Gweon, D. Hansen and E. Perepelitsky at Santa Cruz, and 
with P.W. Anderson, A. Georges, A. C. Hewson, G. Kotliar and Y. Kuramoto elsewhere. This work was supported 



25 



by DOE under Grant No. FG02-06ER46319. 

Appendices 



A. ATOMIC LIMIT t = J = 

A. Double occupancy interpretation of A from the atomic limit 

In order to understand the role of A we study the atomic limit t, J — ^ where this parameter can be introduced 
into the physical Greens function in the form: 

1 — A — 

Q[\ iuj n ] = -2-, (120) 

lUJ n +/X 

and study its dependence on A in the interval [0, 1]. The chemical potential can be calculated from the sum rule 
on the density n of the number of particles N with n = N/N s and temperature T as 

Thermodynamics tells us that the entropy S can be expressed as 



S(n) = -N s I dn' (122) 







and since we know \i from Eq. ( |121[ ) we obtain with y = (1 + A) 

S(n,\) 1 



k B N s 1 + A 



{log4-2/logn-(2-2/) log (2 - y)} . (123) 



we see that its A derivative: k * N |^ = ^ X ) 2 [f + 1°S (1 — f)] * s ne S a tive definite. Thus we see that the entropy 
at a fixed density interpolates monotonically, between the free Fermi limit and the infinite U limits as A ranges 
over its domain < A < 1. The maximum allowed density is reduced from 2 to and thus at A = 1 we have 
a maximum of one electron per site- as expected physically. Thus increasing A from zero effectively removes the 
available states contributing to entropy, its role may be viewed as that of (continuous) removal of states. Thus for 
the equations of motion it is somewhat analogous to the role of Gutzwiller's parameter g in his projection operator 
Yl i [1 — (1 — g)n^nji\ at the wave function level. 

In the atomic limit we can also calculate the entropy at a fixed density of doubly occupied sites d = ^- ^ n it n ii 

as 

^4^ = -d log d-(n- 2d) log(J-d)-(l + d-n) log(l + d-n). (124) 
kb-Ns 2 



An uncorrelated system corresponds to d = where the entropy Eq. (124) is a ma ximu m, while d = for the fully 
projected t-J model. Comparing the two expressions for entropy Eq. ( 123p and Eq. (124), we can express A in terms 



of d at any density. We have thus demonstrated that A is a conjugate variable to the double occupation density in 
this limit. Their explicit relationship is illustrated in Fig. (J6|. 

B. Expansion in A in the Atomic Limit: 



In the atomic limit we set t — >• and J — » so that Eq. ( 65 ) and Eq. ( 66 ) becomes 



Slhf] = go[*,/;Ma] 

= 5[i,/](l-A 7 W)-AM6g[i,f]-Mf,/]. (125) 

Here we split the chemical potential into two pieces fj, = fj, a + A/l^. Thus in this limit g is the free Fermi Greens 
function independent of A, and \i a is the free value i± a — » /xq, the latter determined from the non interacting theory 
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n= .25, .5, .75, 1. 




FIG. 6: The parameter A is determined in terms of the double occupancy d at various densities in the atomic limit as described 
in the text. The arrow indicates increasing density n. Note that the parameter d is scaled into the unit interval. 



in terms of the number of particles. If we turn off the source V the Fourier transforms can be taken as 

r. , r. 1 



■ lUJr, 



H[lUJ n 



l-A? 



lUJ n + Mo 

A \i h g[iu; n ]/i[iu; n ], 



1-Af 



1 + Xfi b g[iu n ] ' 



Thus the physical Greens function 



G[iuj n ] = 



1-Af 



iuj n + Mo + A// 6 

We fix the chemical potentials from the number sum rule as usual and thus 



n 
2 
n 
2 



1 



1 + e~^o ' 



2 7 1 -f e -/5(/x +A/x b ) ' 

We may then solve for /l^'s in terms of the density and obtain 

Ai = k B T log(-^-) 
2 — n 



A Mb = fcsTlog(; 



2 - n 



2- (1 + A)n ; 

Thus the chemical potential Mi has a power series representation 



A" 



m + 1 



2-n 



m+l 



(126) 



(127) 



(128) 



(129) 



(130) 



We see explicitly from Eq. (130) that the A expansion of the the atomic limit is an expansion in An/ (2 — n) i.e. a 
density expansion as well. 



27 



B. THE LOW ORDER CALCULATIONS OF GREENS FUNCTIONS 

A. Greens function to O(A) 

We evaluate the complete starting point of the hierarchy. We start with terms of 0(A°) and end with and 
[g _1 ]i, which are the seeds for the 0(A) terms. 

1. Seed terms and Initialization 



go 1 !*,™] = {[(fi-d Ti - jj )l - Vi]5[i,m] +t[i,m] - V i?m } 

[/i[i,/]] = 1 (131) 



Derived objects 



Jo = S {k) [i,i] 

b[i,m]} = g {k) [m,i] 

te[i,m;i]] = 

[^[i,m;i,fc]] = 0. 



(132) 



[Yi[i,m]] = t[i,m] (g w [i,i] + \e w [m,m])-8[i,m] ^ g (fe) [j,j] - g (fc) [i,j]) 

[*[i,m]] = 0. (133) 

[ A££[i,m;j]] = <W=<W d *M 

[A££[i,m; j,fc]] = <5[*,j] «[m,A]. (134) 

Thus 

[$[i,m]] = 5[i,m] t[i,j] g (fc) [j,i] + i t[i,m] g< fc >[m,m] + ^[i,m] g (fc) [i,j] - ij[i,m] g< fc >[i,m] (135) 

Combining the two we get 

[ Fi [t, m]] + m]] = S[i, m] t[t, j] (g( fc > [t, j] + gW [j, i]) + t[t, m] (g( fc > [i, i] + g( fc > [m, m]) 

-ij[i,m] g< fc >[i,m] - *[i,m] ^J[i,j] g^[j,j] (136) 

Stepping up and final Greens function to O(A). 
To first order in A we collect the above results to obtain the Greens function 

[g- 1 [i,m}] 1 = g(*)[t,»] • V i>ro + <5[i,m] g( fe )[a,i] • vg +<5[i,m] g«[j,j] + ij[i,m] g< fc >[i,m] 

-t[i, m] (g (fc) [i, i] + g (fe) [m, m)) - 5[i, m] t[i,j] (g<*> [i, j] + g<*> [j, i]) (137) 

and the caparison factor: 

UiM\i = -g (fc) [i",i] (138) 
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The FT's of these on turning off the sources is found using g[i , i] — > ^ as 



~ 1 [k}] 1 = ne k + ^J2 Jk -i SM+ \\ J ° n + 2 J2 £ i KM ) 



n 
2 



(139) 



(140) 



Let us note that under the shift Eq. ( jllj ), and reinstating — >• — iaq, the first order correction [g -1 ^]]-,^ shifts 
by 2rm £ + f^j — Sj^o- Therefore this term is invariant under the Shift theorem- II and also the Shift theorem- (LI) as 
specified in Eq. (29). Having verified this, we can drop the terms in the second bracket since it is independent of k. 



B. Greens function to 0(A 2 ) 
1. \i derived objects 

We next start with seed terms of O(A) calculated above and end with [ji\2 and [g _1 ]2- 

Hh m }}i = -g ik) [i~,i] S[i,m] 
Let us calculate the derived quantities from the above at the same level: 

[7[i]]i 



^ (fc) [a,i] 



; (fc) [i,a] = -g[i,i] -g (fe) M 



[Y[i,m]]i 

Zero source Fourier transforms: 



M (fe) [a,i] 



-t[i, m] 



• g (fe) [m,a] 



-g (/c) [ra,z] 



i] • g (/c) [i, i] + -g[m, m] • g (fe) [m, m] 



+ ^[i,m] (j[i,j] g[j,j] -g (fe) [j,j] -f[t,j] g[j,j] -g^J]) 



"T 



3n 2 / n 2 _ n v-^ r i 
-T^Jo+j^e, g[q] 



(141) 



(142) 



(143) 



2. \i derived vertices 



Next we calculate (using lowest order functional derivatives) 



l cr 3 CT4 L ' 

;ZC£[i,m;a] 



iCr2( j3 U/j g,CT4CTl L«^5 « 



-J[z,ra] CTiCT 2 g, 
-S[i,m] aio-2 5 g[i,a] g[a,- t 



At zero sources so with V —> 



' U(a)[i,m;a]] 1 = -2 S[i,m] g[i, a] g[a, i], 
= - <5[z,ra] g[i,a] g[a,i]. 



Z/ s) [z,m;a] 

The zero source Fourier transforms are as follows: 

[K^lPuP^] ] x = - CTKJ2 ^aja^ Yl g ^ gfe + ^ ~ Pi] 

u (a) \pi,P2]\ = -2j2s[q]s[q + P2-Pi] 



(144) 



(145) 



(146) 
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Similarly we find for the four index vertices: 

[ Ksal [h ™>] a , b \ ] i = -6[i,m] (Ti02 &<t 2i <T3 S *A,<Ti gM] g[M], 
U {a) [i,m;a,b] ^ y^o = -2 5[i,m] g[i, a] g[6, i], 

The zero source Fourier transforms are as follows: 

[^3a4 2 [Pl^2;P3,P4] ] x = - 0102 Sa 2ia3 S* 4iai S pi + p4iP2 + p3 g\p 3 ] g\p*\ 

U {a) [p 1 ,p 2 ;P3,P4\ = -2 S Pl + P4iP2 + P3 g[p 3 ] g[p 4 ] 



(147) 



(148) 



5. to O(A) 



We compute from these. 



[*(*01i = ^\£p + \ £ k + \ J k-p\ g[p] ^(p,*)]^^^-* g M [ w(o) (p» fc ;9+P-fc,9) 



= - ^ (e P + £k+ q -p + £fe + £ q + Jfe-p) g[p] g[g] gfe + fc - p] 

p,q 

= ^2W(k,q;q + k-p,p) g[p] g[q] g[q + k - p] 



(149) 



4- Stepping up: \i to 0(X 2 ) 

Stepping up, we calculate 

WM\2 = S[i,m] g[r,z]g (/c) [r,z]+ [^m]^ 

Hence at zero sources, the Fourier transform reads 

^2 



WW2 = ~i 



+ £k+ q -p + e k + e q + J k - P ) g[p] gM g[« + & ~ p] 



(150) 



(151) 



Note that \ji\2-, are invariant under all three shift theorems. It is clear that this is a more non trivial application 
of the theorems than those in the lowest order. 



5. g 1 derived objects 



Let us now start with: [g 1 [i,m]] 1 given in Eq. (137): 

[g" 1 [i, m}} 1 = [g<*> [i, i] ■ V iim + S[i, m] g( fe ) [a, i] ■ V$ 

-t[i, to] (g( fe > [i, i] + g (fe ) [to, to]) - 5[i, to] t[i, j] (g (fc) [i, j] + g (fc) [j, 



S[i,m] ij[*J] g (fc) [j,j] + ^J[*,m] g (fe) [i,m] 



(152) 



Three point vertex: 



6. Vertex functions to O(X) 



[A--[i,m;a]] 1 = -(^^) [g^M^ 



[(/) + (//) + (!!/)]£ 



C^2 
Cr3(T4 5 



(153) 
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where the terms (J), (II), (III) refer to the three square bracketed terms in Eq. (152). 
For the first term we calculate 



( T\CTlCT2 

K 1 Jo 



v^o = 0. 

For all other terms we can use a simple calculation: 



(154) 
(155) 



Therefore 



(criCT2)ga 2 a 3 [p, rjg^ [r, g] 
Af (a) [p,g;r] v^o = 2 g[p,r]g[r, g] 



(7/) (a) [z,ra;a] = -£[z,ra] J[z,j] g[j,a] g[a,j] - J[z,ra] g[z,a] g[a,m], 



(156) 
(157) 



and 



(m) (a) [z,ra;a] = 2[t[z,ra] (g[z, a]g[a, z] + g[m, a]g[a, m]) + £[z, m] £[z, j] (g[z, a]g[a, j] + g[j, a]g[a, i]) ] (158) 
Zero source Fourier transforms read as: 

(II) {a) [p 1 ,p 2 \ = -J P2 - Pl ^2g[q] g[q+P2 -Pi] - J Pl -qg[q]g[q + P2-pi] 



(III) {a) [p 1 ,p 2 \ = -2^2g[q]g[q+p 2 -pi] {e P2 + e Pl + e q+P2 - Pl + e q } 



(159) 



Hence adding up we obtain: 

[ KlallPuP^] ]l = -(^2)Sa 2 a 3 Sa 1<T4 ^ *pi+p 4 ,P2+P3 g^g^] + £p 2 + £p 3 + £ P4 + \ ( J Pi~P2 + J Pi-p 3 )} 



^l^Wa W 4 gfe]g[P4] [W(p 2 ,P3,P4,Pl) + W(p 2 ,P3,Pl,P4)] 



P3,P4 



A (a) [pi,p 2 ] 



5^ g[P3]g[P4] [W(p 2 ,P3,P4,Pl) + ^(P2,P3,Pl,P4)] 



P3,P4 



Note that rotation invariance relations imply that since [A^fpi,^] ] x = 0, we must have 

AWbl.P2]l = -W t) \Pl,P2}} =1 W a) \pi,P2] 



1 2 



(160) 
(161) 



(162) 



The four point vertex. The calculation proceeds similarly: 



[A^[i,m;a,6]] 1 = -(™^) [g^ 2 M]i 



[(JV) + (V) + (VJ)]££ 



(163) 



Here the terms (/V) — (V/) refer to the three square bracketed terms in Eq. (152). For the first term we calculate 



a.b 



= -5[i, a]5[m, b] [(<7i<73)<W 4 8*3*1 & 1 ~ a M^> m ] K^i^) ^ 2 ,a 3 g^^ [&, «]] 
(^)a 3 a4 v^o -5[i,a]5[m,b] S ai(T3 J a2a4 g[z,z] - 5[i, a] 5[i, m] (a^ 2 ) S ai&4 S a2&3 g[b,i] 



(IV) 



(a) 



6[i, a]5[m, b] g[z, z] — 2 5[z, a]5[z, m] g[b, z] . 



(164) 



(165) 
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This term is seen result in a violation of Eq. (47) and Eq. (70) for reasons discussed there and in the second remark 



below Eq. (66), and therefore is dropped below. We have carried it in the calculation, and demarcated it with the 



underbrace, in order to see its (minor) contribution explicitly before dropping it. 
For all other terms we can use a simple calculation: 

yilZ [p, q; r, s] = (^A^)gW a [p, q] 
y ia) {p,q;r,s] v ^o = 2g[p,r] g[s,g] 

Therefore 

(V") (a) [z',ra;a,6] = -5[i, m] J[i, j] g[j, a] g[6, j] - J[z, m] g[i, a] g[6, m] 



(166) 



(167) 



(^/)W[i,m;a,6] = 2[£[z,m] (g[z, a]g[6, z] + g[m, a]g[6, m]) + S[i, m] t[i, j] (g[i, a]g[6, j] + g[j,a]g[M])] (168) 
(jy)( a ) [i 5 m; a, 6] = a]5[ra, 6] g[i, z] — 2 a]5[i, m] g[b, z] 
(V) (a) [z,m;a,6] = -6[i,m] J[z,j] g[j,a] g[6,j] - J[z,ra] g[z,a] g[6,m] (169) 

(V) (a) [pi,p 2 ,P3,P4] = - (Jp 2 -pi + Jpi-ps) g[P3] g[P4] 

(l//) (a) [z, m; a, 6] = 2 [ t[z, m] (g[z, a]g[6, z] + g[m, a]g[6, m]) + £[z, m] t[z, j] (g[z, a]g[6, j] + g[j, a]g[6, z]) ] 
(y/) (a) [pi,p 2 ,P3,P4] = -2{e P2 + +£ P3 + £ P4 } g[p 3 ] g[p 4 ] S Pl + P 4 iP2 + P3 (170) 

Hence 

~ ^Pl,P3^P2,P4^O-i<73^<T2O- 4 g[0 ] (<7l<72) ^cria 4 ^cr 2 a3 ^Pl+P4,P2+P3 (l^l) 



Thus 



A (a) [pi,p 2 ,p 3 ,P4] 



— 2 ^pi+p 4 ,p 2 4-p 3 g[p3]g[P4] {^pi + + £ P3 + £ P4 + 2 (^Pl _ ^2 + ^Pl-P3)} 

+ ^Pl,P3 <Wp4 g[0~] _ 2 ^Pl+P4,P2+P3 g[P4] 



(172) 



Comparing Eq. (160) and Eq. (172), we see that other than the term with underbraces, these vertices satisfy Eq. (47) 
or Eq. (PrO), 



7. $ to O(A) 



We now assemble terms: 



P \ " / pgr 



(173) 
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Let us rewrite this as (k — > P2,p — > Pi, q — > Pi) 

Pl ^ ' 

= T^ £ P3g[P3] - epi+p4-p 2 g[Pl]g[P4] 



E o £ ps sbi] A (a) bi 



i ^ 2 

Pl+P4=P2+P3 



0l,P4 



"2 E gbl]gb3]g[P4] U P1 + + \j Pl - P2 j + £p 3 +^4 + ^ (Jpi-P2 + J P1-P 3 )I 

E g[Pi]g[P3]g[P4] £ P3 {e Pl S P2 £ £>3 "I" £ P4 "I" o (Jpi~P2 Jpi-Pz)} (1^4) 



0l+P4=P2+P3 



The first line with underbraces arises from the term A in Eq. (|169|), or Eq. (11651) and Eq. (172) which disobey the 



relation Eq. (47) or Eq. (70). It gives a static but momentum dependent contribution, and we will drop it as discussed 
below Eq. (47) and in the second remark below Eq. (66). The rest are combined and rearranged to give 

[*(P2)]i = - E g[Pi]g[P3]g[P4] (e 

Pl ^P2 £ P3 ^P4 Jpi—P?) {^Pl £ P2 £ P3 ^PA o (Jpi—P2 Jpi~ P3)}' 

Pl+P4=P2+P3 

= \ E %\Pi\%\Ps\z\pd w fa^ ( 175 ) 



Pl,P3,P4 



where in the first line we symmetrized further in pi O p±. 

We can bring this into standard notation by sending P2 —> k,pi p^ps — > q^p^ —> k + q — p: 

M k )\i = ~E g M g bl g[& + 4-p] 

X (gjfe + £p + £ g + Ek+q-p + {^/c + £p + Sq + Ek+q-p + - (Jfc-p + ^p-g)} 

= ^Z)sM gW g[* + 9-p]W^(fc,9;9 + fc-p,p)[W^(fc,9;9 + fc-p,p) + W(fc,g;p,g + £;-p)] (176) 



9. Stepping up and final Greens function to 0(A 2 ). 



We are now in a position to put together the second order result for g _1 and also \i. Recall that [g -1 [&]] 2 = 
— [Yi[k] -\-<$>[k]] 1: where these variables are calculated in Eq. (143) and Eq. (176). Hence we can now compile the 
equations of the second order theory with sources turned off: 

g" 1 ^) = go 1 (A : ) + A[g- 1 (fc)] 1 + A 2 [g- 1 (fc)] 2 + 0(A 3 ), 



[g"M% 
[g" 1 M] 2 



iuj n + M - e k - 1 J 



n s k 



\ Jk -* g M + ( \ J ° n + 2 £ <i g M ) 
q V q J 

Tlj^ 1 ( 77/ 2 77/ 1 

— s k - Mk)}, - —Jo + - Y, £ i sk] 



(177) 



We thus see that all the computed [g x )j are invariant under the three shift theorems. Having verified this, we can 
drop the constant term and adding up terms to 0(A 2 ): 



/ 3r? 2 \ 1 

g" 1 ^) = iw n + M - f 1 - A n + A 2 —j e k + A ^ ^k-q gfe] - A 2 + 0(A 3 ), 



(178) 



along with the definition in Eq. (176). 
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Similarly the caparison factor /i is found to 0(A 2 ) as: 



fi[k] = 1 + A [nlk^ + X 2 [fi[k]] 2 + 0(A 3 ), 



n 
2' 



(179) 



Adding up terms to 0(A 2 ) we obtain: 



H[k] = 1-A| + A 2 ^ + A 2 [*(fc)]i + 0(A 3 ). 



(180) 



along with the definition in Eq. (149). 
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